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Disorder modifies the sound-wave excitation spectrum of Bose-Einstein condensates. We con- 
sider the classical hydrodynamic limit, where the disorder correlation length is much longer than 
the condensate healing length. By perturbation theory, we compute the phonon lifetime and the 
correction to the speed of sound. This correction is found to be negative in all dimensions, with 
universal asymptotics for smooth correlations. Considering in detail optical speckle potentials, we 
find a quite rich intermediate structure. This has consequences for the average density of states, 
particularly in one dimension, where we find a "boson dip" next to a sharp "boson peak" as function 
of frequency. In one dimension, our prediction is verified in detail by a numerical integration of the 
Gross-Pitaevskii equation. 
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I. INTRODUCTION 

Disorder is a key feature for the understanding of the 
properties of matter. In a disordered environment, waves 
can become coherently localized, leading to suppression 
of transport [1, 2]. Also strong interaction can produce 
an insulator, via the Mott-Hubbard transition [■'>]. The 
combined effects of interaction and disorder, despite be- 
ing studied for decades, still hold surprises. Here we are 
interested in the influence of a spatially correlated disor- 
der potential on the low-energy excitations of an interact- 
ing Bose-Einstein condensate (BEC). As Goldstone exci- 
tations, these low-energy Bogoliubov excitations feature 
a linear, phonon-like dispersion relation LUk = ck with 
sound velocity c. The sound velocity is of particular in- 
terest because it determines the range of superfluidity, 
according to the Landau criterion, and it determines the 
density of states, which enters virtually all physically rel- 
evant quantities. Moreover, the speed of sound is directly 
measurable in cold-atom BECs [ ], where well-controlled 
optical speckle potentials with interesting spatial corre- 
lations can be studied [2, ■)]. 

Calculating the effective speed of sound c in disordered 
systems is far from trivial, with different approaches lead- 
ing to different predictions. Perturbation theory, on the 
one hand, predicts an increased speed of sound due to 3D 
uncorrelated disorder [(i, 7]. On the other hand, within 
a self-consistent nonperturbative approach Yukalov and 
Graham [-'] reported numerically a slight decrease. For 
disordered hard-core bosons on a lattice, Zhang [ ] found 
a decrease of c to fourth order in disorder strength, with- 
out information on the second-order effect. 

To clarify the situation with suitable parameters for 
present-day BEC experiments, we study in this article 
phonon excitations of a BEC in the strongly interact- 
ing case where the chemical potential fi is much larger 
than the disorder strength V, and where the conden- 
sate healing length ^ — h/ ^1ra[i is much smaller than 
the disorder correlation length a (see Fig. 1). Without 
disorder, the BEC is contained in a very shallow trap 
and has a constant density no = [ijg in the region of 
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FIG. 1: (Color online) Schematic Id representation of the 
system under study; an interacting Bose-Einstein condensate 
with original homogeneous density no = ^/g (dashed black 
line) is exposed to a weak, spatially correlated random po- 
tential V{r) (solid blue), here a blue-detuned speckle poten- 
tial with amplitude V = O.I/.1, centered on the mean V = G 
(dotted blue). We consider the Thomas- Fermi regime where 
the healing length ^ is much shorter than the disorder cor- 
relation length a. The resulting ground-state density (solid 
black) [eq. (9)] mirrors the disorder while leaving the total av- 
erage density and number of particles constant. On top of this 
disorder-modified ground state, an elementary plane- wave ex- 
citation (green, plotted around 1) propagates with wave vec- 
tor k, here with ka = 1. We calculate its effective speed of 
sound and the corresponding average density of states. 

interest. In the presence of smooth disorder with corre- 
lation length cr ^ ^, the BEC ground state density fol- 
lows the external potential with the Thomas- Fermi pro- 
file no(r) = no[\ — V{r)/n]. A long- wavelength density 
deviation 5n(r,t) from this ground state obeys the wave 
equation 

\y ■c^{r)V -dl]5n^Q (I) 
1/2 

where c(r) = c[\ — V{r) / n] ' is the local speed of 
sound, deviating from the clean value c — \J ^/m = 

■y/ gno/m. This is a prototypical wave equation in a 
medium with random elasticity, but constant mass den- 



2 



sity [10, 11]. Quite often, the opposite case is studied, 
with random masses and constant elasticity, or equiva- 
lently, a fluctuating index of refraction [12, 13, 14]. The 
disorder potential may always be taken at zero average 
V{r) — 0. Its strength is characterized by the variance 
V{r)^ — V'^, and we suppose weak disorder with V <C /i. 

Consider now a sound wave with wave vector k evolv- 
ing on the disordered potential background with corre- 
lation length (7 (see Fig. 1). If the wavelength is much 
longer than the correlation length, ka <C 1, the excita- 
tion averages over the potential fluctuations and, to a 
first approximation, it seems reasonable to replace (1) by 
its ensemble-average [12]. But then we have no net ef- 
fect on the speed of sound since c'^{r) — (? , exactly. If, 
on the other hand, the wavelength is much shorter than 
the correlation length, fccr ^ 1, the excitation evolves in a 
locally constant potential, which should result in an aver- 
age speed of sound c = c[l - V{r)/^Y/'^ « c[l-|y^//i^]. 
It turns out, however, that both these naive reasonings 
fall short. 

In order to give the correct answer right away, our main 
results are briefly summarized in the following section 
II. Section III then presents the general hydrodynamical 
perturbation theory, from which detailed results on the 
speed of sound are derived in section IV. In section V, we 
analyze the implications of these results for the disorder- 
averaged density of states. A short conclusion together 
with a brief comparison to related works are contained 
in section VI. 



II. MAIN RESULTS 

The effective speed of sound in a disordered interact- 
ing Bose gas, properly defined as oJk/k = c from the 
single-excitation dispersion relation, is affected by scat- 
tering processes via virtual intermediate states such that 
a purely local description fails. We find that the correc- 
tion Ac — c — c to the speed of sound has in d dimensions 
the limiting behavior 
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ka > 1. 
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These limits imply that the curves for different dimen- 
sions have to intersect around fccr ~ 1 (see also Fig. 2 
below). The precise shape of Ac/c at intermediate ka 
depends on the details of the disorder pair correlation 
function. But clearly, there is a negative correction, of 
order F^Z/i^, in all dimensions and for any disorder with 
finite correlation length cr 3> ^. 

A reduced speed of sound implies that the free density 
of states (DOS) of single excitations. 
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is replaced by an enhanced average density of states (AV- 
DOS) 'p{^)- Our results for this disorder-induced correc- 



tion can be cast into the form of a function 

gdiyjojc) = \p(uj) - po(^)] /Po(t^) (5) 
that depends only on the reduced momentum n = ua /c: 
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K < 1, 
K > 1. 

(6) 

Gurarie and Altland [12] suggested that one should be 
able to deduce from the asymptotic values and the curva- 
tures of such a scaling function whether the AVDOS ex- 
hibits a "boson peak" at intermediate frequency to « c/a. 
The asymptotics of the scaling function in our case allow 
for a smooth, monotonic transition between the limiting 
values in any dimension d. Thus one has no reason to ex- 
pect any extrema in-between, which is indeed found to be 
the case in two and three dimensions. In d = 1 however, 
we find, by analytical calculation for the experimentally 
relevant case of an optical speckle potential, a quite non- 
monotonic AVDOS with an intermediate dip followed by 
a sharp peak at ua/c = 1. 

III. CLASSICAL HYDRODYNAMIC THEORY 

We start our detailed analysis of the mean-field BEG 



order parameter ^I* 



/n e 



in terms of the hydrody- 



namic variables condensate density n = |^'|2 and phase 
(/), which determines the superfiuid velocity v = ;^V0 
[15, IG]. The grand-canonical Gross-Pitaevskii energy 
functional for the BEG in presence of an external poten- 
tial V{r) is 

2m 



d'^r 



(V^^)%7^(V(/))^ 



(7) 



+ {V{r)-fi)n+^n' 



The saddle-point equations SE/dn\o — and SE/64)\o ~ 
imply that the ground state has constant phase 0o or 
zero superfiuid velocity Vq — and a density no(r) that 
obeys the stationary Gross-Pitaevskii equation 



^vVMr) 



2m 



+ gnQ{r) = fi-V{r). 



(8) 



We now restrict our analysis to the case where the heal- 
ing length ^ is much shorter than the disorder correlation 
length a. In this regime, the quantum pressure in (8) is 
negligible, and the external disorder potential leaves the 
Thomas-Fermi imprint 



(9) 



This solution is also directly obtained if one drops the 
density-gradient contribution in (7) to use 

|^n(V^)2 + [y(r)-~/i]n+ V|. 

(10) 
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Formally, this formulation corresponds to the limit ^ — > 
0, and all further results can only depend on the reduced 
momentum ka [17]. 

The speed of sound characterizes the dynamics of small 
deviations 6n{r,t) = n{r,t) — no{r) and 6(j){r,t) — 
<j){r, t) — (j)Q from the ground state in the long- wavelength 
regime fc^ <C 1. We can therefore develop the energy 
functional (10) to second order around the ground-state 
solution, E' = E'q + F'[Sn,S(f)], to obtain the relevant 
quadratic energy functional 



F'[Sn,6(j)] = 



^ / d^r|^no(r)(V<50) 



2 



gSn' 



(11) 



Importantly, the external disorder potential has shifted 
the ground-state solution according to (9), around which 
we now consider the dynamics of fluctuations. Density 
and phase are conjugate variables with the equations of 
motion 



hdtSn 



SF' 



'hdtS<f> - 



SF' 



(12) 



In terms of density and superfluid velocity, they read 

dt5n + V ■[nQ{r)v] = Q, (13) 

dtv = -—V6n, (14) 
m 

and are recognized as the linearized versions of continuity 
equation and Euler's equation for an ideal compressible 
fluid, respectively. These can be combined to a single 
classical wave equation 



(15) 



This equation is equivalent to (1), but now written in a 
form amenable to systematic perturbation theory for a 
weak external disorder potential V{r). 



A. Perturbation theory 

Translation invariance of the free equation suggests to 
use a Fourier representation in space and time. 



The single-excitation dispersion relation can be derived 
from the corresponding Green function. The free Green 
function is diagonal in fc. 



(18) 



Taking the disorder average of the full Green function 
G — [Gq^ — V] leads in the standard way to [14, 19] 



G{k,uj)^ [Go{k,Lj)-' ~nk,io)] 



(19) 



The poles of this average Green function at = c^k"^ + 
Y.(k,uj) now determine the effective dispersion relation. 
The so-called self energy S(fc, w) is given to leading order 
in disorder strength by the Born approximation: 



E(fc,tj) 



d'^fc' 



[VV]fc^,,Go(fc',^). (20) 



The scattering potential correlator 

Wkk' = [fe • fc'] ' <^''Pd{\k' - k\a) (21) 

involves the dimensionless fc-space correlator of the bare 
potential 



(22) 



Its real-space correlator Cd{r/(j) — V{r)V{Q) jV"^ is as- 
sumed to be isotropic. We will consider correlated poten- 
tials for which G{r/a) decays from Cd(0) = 1 to on the 
length scale a. The smoothness of V{r) implies that the 
power spectrum Pdin) decreases rapidly to as function 
of K = ka. 

Applying Sokhotsky's formula {x + iO)^^ = P - — 
iTrS{x) to the free Green function in (20), we can evaluate 
the real part and the imaginary part of the self-energy 
separately. The imaginary part determines the lifetime 
of the excitations, whereas its real part shifts the 
speed of sound by Ac: 



T,{k,ck) _ Ac 
2c2fc2 ~ ~ 
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(23) 



To leading order in V, the on-shell dispersion ui = ck is 
used for evaluating the self-energy. 



/d'^k' 
j^-^^kk'Srifs' . (16) B. Scattering rate and Id localization length 



The disorder potential causes scattering k ^ k' oi plane 
waves with an amplitude 



V, 



fcfe' 



— {k-k')Vk-k'- 



(17) 



The factor k k' originates from the mixed gradient in (15) 
and implies pure p-wave scattering of sound waves [18] 
in contrast to s-wave scattering of independent particles 
[19]. 



Calculating the imaginary part in (20), the scattering 
rate at frequency u — ck can be expressed as 

7(^) = ^u:^P^{LoVfd {ojajc) . (24) 

The last factor is the angular average of the correlation 
function on the energy shell: 

/d(K) = y df^d(cos^)2Pd(2Ksin^). (25) 
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The squared cosine under the integral goes back to the 
[k ■ fe'] in the potential correlator (21), being again 
characteristic for p-wave scattering of sound waves. 

In one dimension, there are only the two contributions 
6 — 0,7r of forward- and backscattering, respectively, 
such that 

li^) - [Pi (0) + Pi (2^a/c)] . (26) 

We note in passing that the Id backscattering process 
k —k described by the second contribution Pi (2k) is 
known to induce strong, Anderson localization of the ex- 
citation in the disordered potential [20]. The backscatter- 
ing rate is directly proportional to the inverse localization 
length Floe = 7bs/2co describing exponential localization 
[21]. Taking the backscattering contribution of (26), we 
find 



Pi(2fca), 



(27) 



which agrees with the findings of a hydrodynamic the- 
ory similar to ours [ ] and also with the sound-wave 
limit of Bogoliubov excitations considered in [ ] . It 
should be noted that these latter approaches employ the 
phase-formalism that is particularly suited for Id sys- 
tems, whereas our Green-function theory permits to go 
to higher dimensions without conceptual difficulties. 

In any dimension, the phase function fd^K) in (25) 
tends to a constant for small k, and the scattering rate 
of low-energy excitations rate tends to zero as 7/w oc w''. 
Also the localization length in Id diverges as F^"^ oc uj~^ 
at low frequency. In higher dimensions, it is known to 
be even larger, if not infinite [L!]. This assures that low- 
energy excitations are long-lived and extended. It is thus 
meaningful to discuss their effective sound velocity. 



IV. EFFECTIVE SPEED OF SOUND 

From (20) and (23), the speed-of-sound shift Ac in 
any dimension d is obtained as a Cauchy principal value 
integral over the potential correlation: 

Ac ^ IV\ f d'^ik'T) [k ■ k'] ' Pd{\k' ~ k\a) 
c 2 J {2nY k^ (fc'2 - fc2) ■ ^ ' 

In the limits fcu 3> 1 and fccr <C 1 where the potential 
appears very smooth or (5-correlated, respectively, over a 
wavelength of the propagating excitation, this correction 
is independent of the precise form of the bare potential 
correlator Pd{K) (see (2) and (3) above and the detailed 
derivation in section IV C below) . Let us then discuss the 
interesting, detailed form of this correction as function of 
K — ka in d = 1. For concreteness, we study the case 
of an optical speckle potential, which has recently been 
successfully used in experiments on Anderson localization 
of matter waves [ ']. 



A. Speckle potential 

By focusing a laser beam through a diffusor, the con- 
densate is subject to a random lightshift potential pro- 
portional to the intensity of the laser field [ J4] . The one- 
point potential value V{r) of a speckle pattern [25] has 
the skewed probability distribution 

P{w)dw ^ e{l + w)exp[~{l + w)]dw, (29) 

for w = V{r)/V. For this one-sided exponential, odd 
moments, such as V{r)^, are different from zero. A blue- 
detuned light-shift potential with y > features repul- 
sive peaks (this case is depicted in Fig. 1), whereas a 
red-detuned one with V < consists of attractive wells. 
As far as spatial correlations are concerned, the laws of 
optics forbid variations on a length scale shorter than the 
correlation length a, which depends on the laser wave- 
length and the geometry of the imaging system, but typ- 
ically ranges around 1 /^m. In one dimension, the corre- 
lation function is 



Pi(K) = -(2-|K|)e(2-|K|). 



(30) 



Its bounded support in fc-space implies that within the 
Born approximation, backscattering and inverse localiza- 
tion length vanish for ka > 1; however, exponential local- 
ization still prevails due to higher orders in perturbation 
theory ['>]. 



B. In dimension d — 1 

The principal- value integral (28) over the piecewise lin- 
ear function (30) is elementary, and we find a speed-of- 
sound correction at k = ka of 
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(31) 



Its limiting values are Ac/c = —^V'^/fi'^ for small k, 
and Ac/c = —'^V^/fj,^ for large k, as stated in (2) and 
(3). This correction to the speed of sound, plotted in 
Fig. 2 as function of k = ka, shows a rather intricate, 
non-monotonic behavior. Notably, there is a logarith- 
mic non-analyticity at k = 1, the value beyond which 
backscattering is suppressed. The speed-of-sound correc- 
tion is clearly negative for all k, which may come as a 
surprise in view of [6, 7] . 

In order to check this prediction in detail, we have 
numerically integrated the full Gross-Pitaevskii equation 
describing an elementary excitation with fixed k on top 
of the numerically determined groundstate in a speckle 
potential with V — ±0.03/i and variable a. This simu- 
lation operates at a small but finite value of k£_ ~ 0.05 
and includes the full quantum pressure. Moreover, it 
does not rely on a linearization for small excitations nor 
perturbation theory in V. We extract the effective dis- 
persion LOk by monitoring the phase of ^'fc(t) and then 
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FIG. 2: (Color online) Relative correction to the speed of 
sound as function of reduced momentum. Black line: analyt- 
ical prediction (31) for d = 1 in the limit ^/ct = 0. Blue/red 
symbols: data from a numerical integration of the Gross- 
Pitaevskii equation (with small but finite healing length, such 
that A;^ = 0.05) at ?; = V/ ^ = ±0.03, respectively, averaged 
over 50 realizations of disorder. In the shaded area, the data 
strays far from the analytical prediction because its condition 
of applicability cr ^ ^ becomes invalid. Inset: Theory (28) 
for d = 1 (black), d = 2 (violet), and d = 3 (orange) together 
with the limiting values (2) and (3). 



find c = Wfe/fc. As shown in Fig. 2, the data agree 
beautifully with (31) in its realm of validity, cr 3> ^. 
When the correlation length decreases towards the heal- 
ing length (shaded area in Fig. 2), the correction van- 
ishes at fixed disorder strength, because the condensate 
density is smoothed with respect to the Thomas-Fermi 
profile [-"()]. But in any case, only negative corrections 
are found in d = 1. 



FIG. 3: (Color online) Relative correction to the speed of 
sound at ka = 1, divided by as function of disorder 
strength v = V/^j,, taken from the numerical integration with 
fc^ — 0.05. Dashed line: analytical prediction Ac/(cw^) = 
iln2- I -0.3267. 



nearly everywhere under the integral. Then 



Arf(oo) 



{2n) 
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Sd 
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(33) 



and with / f^Pdiq) = Cd{0) 



the limit ka 



1, we arrive at (2). In 
oo we proceed similarly with 77 ^ 0. The 



angular integrand reduces to 1 -I- [2r] cos P + rj'^] , whose 
principle- value integral evaluates after some algebra to 

Sd d + 2 



Ad{0) 



{2ny 4 



(34) 



which leads to (3). 



C. In higher dimensions 

In higher dimensions, the integral (28) over the corre- 
lation functions (for speckle, see [19]) is sufficiently com- 
plicated that analytical solutions like (31) are not avail- 
able in general. (As an exception to this rule, we find 
for the 2D speckle correlation Ac/c = ~ ^^^) 

for K > 1.) But in all cases, the principle- value integral 
(28) can be evaluated numerically. The inset of Fig. 2 
shows the corresponding curves. Short-range correlated 
potentials {ka <C 1) affect low dimensions more than high 
dimensions and vice versa. 

The limits (2) and (3) can be calculated analytically as 
follows. It is useful to rewrite (28) in terms oi i] = l/ka 
as 

Ac^ f c¥'q Pd{q)[l + vqcosPf 

c 2^2 J (27r)^ 2?7gcos/3-hr?2g2 • ^ > 

Denoting the angular part of the integral by Ad{'>]q), one 
arrives at the radial integral dq q'^~^ Pd{q) Adirjq) . In 
the limit fccr <C 1, the parameter rjq tends to infinity 



D. Numerical investigation beyond Born 

The perturbative analysis relies on the Born approxi- 
mation (20), so that good agreement with the true values 
is only expected at rather small disorder. Could a larger 
disorder strength reverse the sign of the correction? We 
have numerically investigated different values of v = V/ fi 
at fixed ka = 1. In Fig. 3, we show the data divided by 

such that the Born approximation shows as a horizon- 
tal line. As expected, for small \v\ the agreement is very 
satisfactory. One can distinguish a third-order correction 
0{v^) as a linear trend with negative slope; if needed, it 
could be calculated pushing (20) beyond the Born ap- 
proximation [ ]. In a Gaussian model with a symmet- 
ric probability distribution [ ], such a third-order term 
would be absent. 

The error bars in Figs. (2) and (3) indicate the es- 
timated error of the mean after ensemble-averaging over 
50 realizations of disorder. Fig. 4 displays exemplary his- 
tograms of the values obtained for different disorder real- 
izations. Clearly, the probability distributions are single- 
peaked with well-defined averages on the negative side. 
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FIG. 4: (Color online) Histograms over 50 realizations of dis- 
order for the relative speed-of-sound correction for different 
values of 1) = V/^ obtained by numerical integration of the 
Gross-Pitaevskii equation with an excitation at fca = 1. Ver- 
tical lines indicate the average values plotted in Fig. 3. 



That the speed of sound has self-averaging character was 
to be expected since a plane wave samples different spa- 
tial regions at once. At strong disorder with v > 0.1, 
the speckle disorder with its unbounded probability dis- 
tribution is likely to fragment the condensate, and the 
concept of a unique, well-defined speed of sound becomes 
questionable. As a precursor, we already observe a slight 
broadening of the probability distribution for v = +0.1 
(lower right panel). 

From the data shown, we conclude that the correction 
to the speed of sound remains negative over the entire 
interval of interest. 



V. DENSITY OF STATES 

Knowing the speed of sound, we can compute the av- 
erage density of states (AVDOS) 



(27r)° 



(35) 



using the effective dispersion uj^ = c{k)k in the pertur- 
bative limit where j/uj <^ 1. Denoting, similarly to [12], 



p(w) = poiuj) [1 + gd{uja/c)] 
we find for the relative correction 



5d(«) = - 



dn 



Ac 
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(36) 
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with the limiting values (6) [2n]. The scaling functions 
gdi^) are plotted in Fig. 5 for dimension d = 1,2,3. In 
one dimension, 



g,{n)^-[l + -ln 



K- 1 



3k^ 
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FIG. 5: Correction to the density of states (/d(K) = {p~ Po) / Po 
divided by the squared disorder strength v = V/n as function 
of reduced momentum k — coa/c in dimension d = 1, 2, 3. At 
K = 1, the momentum beyond which elastic backscattering 
becomes impossible in the Born approximation, there is a 
logarithmic divergence in d = 1, a kink in d = 2, and a 
curvature discontinuity in d = 3. 



shows a pronounced dip around k « 0.7 and a sharp log- 
arithmic divergence at n = 1. This particular structure 
is a consequence of the Born approximation, more specif- 
ically the non-analyticity of the speckle pair correlation 
function at the boundary of its support. But a local max- 
imum is also found for other correlated potentials with 
fast enough decay of Pd{q) such as Gaussian correlation 
[ ] . The existence of this "boson moat" could not be in- 
ferred from the asymptotics of gi{K) alone [12]. Indeed, 
expanding the asymptotic behavior as 



9d{K) 
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V X 



+ ...), 



+ ...), 



we find = —1 — ijliiK 



< and 



K < 1, 
K > 1, 



(39) 



opposite sign. Together with the fact that is larger 
than (3^ , these asymptotics would be compatible with a 
monotonic behavior and thus are not sufficient to infer 
the existence of intermediate extrema. 

In two dimensions, the scaling function is exactly con- 
stant for K > 1 and thus aj" = which seems to happen 
also in other cases [12]. At k = 1, there is a kink, but 
overall, 52('«) shows a monotonic behavior without local 
extrema. In three dimensions, the logarithmic singular- 
ity has moved to the second derivative of g3{K), which is 
hardly resolvable in the figure, and < as expected 
[12], leaving an all but structureless AVDOS. 



VI. CONCLUSIONS 

In conclusion, we have perturbatively calculated the 
infiuence of a weak spatially correlated disorder poten- 
tial on the sound-wave spectrum of Bose-Einstein con- 
densates in the hydrodynamic limit ^ <g; a,l/k and ar- 
bitrary dimension. The sound-wave lifetimes are long 
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enough to observe a disorder-induced correction to the 
speed of sound, which is found to be reduced. For the 
experimentally relevant case of an optical speckle poten- 
tial, we compute the correction to the speed of sound an- 
alytically. A numerical integration of the full mean-field 
dynamics in d = 1 confirms our prediction in its range of 
validity and even allows to access non-perturbative dis- 
order strengths. 

The present hydrodynamic theory compares well to re- 
sults in d = 1 from the phase-formalism approaches of 
Bilas and PavlofF [2"J] and Lugan et al. We find per- 
fect agreement concerning the localization length, which 
we obtain from the backscattering rate. 

However, our results are in contrast to the impact 
of uncorrelated disorder in three dimensions, for which 
Giorgini, Pitaevskii and Stringari [ ] have predicted a 
positive correction to the speed of sound. Yet, at present 
there appears no contradiction between their results and 
ours. We have used a simple hydrodynamic description 
valid for f ^ cr, 1/fc that cannot cover the case of a 
truly ^-correlated disorder, spatially varying on a scale 
(T ^ ^, considered by Giorgini et al. In particular, for 
such rapidly varying potentials the Thomas-Fermi ap- 
proximation Eq. (9) for the ground state density does 
not hold any more and should be replaced by the solu- 



tion of (8), which then shows the smoothed imprint of 
the disorder potential [■2()]. 

Furthermore, we have determined the average sound- 
wave density of states. In low dimensions, its structure is 
very rich, including a "boson moat" consisting of a broad 
dip followed by a sharp peak in d = 1. As a rule, specific 
correlation-related features tend to be washed out by in- 
tegration in higher dimensional fe-space. Thus we expect 
arguments on general grounds [12] to hold more reliably 
in higher dimensions. Conversely, the low-dimensional 
behavior may escape a bird's-eye view and require de- 
tailed calculations. We have presented such a calculation 
for spatially correlated speckle disorder, so that our re- 
sults should be of immediate use for cold-atom experi- 
ments. 
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